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Abstract. In this paper we present some of the recent progresses in multiple zeta values (MZVs). 
We review the double shuffle relations for convergent MZVs and summarize generalizations of 
the sum formula and the decomposition formula of Euler for MZVs. We then discuss how to 
apply methods borrowed from renormalization in quantum field theory and from pseudodifferential 
calculus to partially extend the double shuffle relations to divergent MZVs. 



The purpose of this paper is to give a survey of recent developments in multiple zeta values 
(MZVs). We emphasize on the double shuffle relations which underlie the algebraic relations 
among the convergent MZVs, and on renormalization methods that aim to extend the double 
shuffle relations to MZVs outside of the convergent range of the nested sums defining MZVs. We 
also provide background on double shuffle relations and renormalization, as well as the closely 
related Rota-Baxter algebras and some analytic tools in pseudodifferential calculus in view of 
renormalization. 

1.1. Double shuffle relations and Euler 's formulas. A multiple zeta value (MZV) is the spe- 
cial value of the complex valued function 



at positive integers si, ■ • • , with s\ > 2 to insure the convergence of the nested sum. MZVs are 
natural generalizations of the Riemann zeta values g(s) to multiple variables. The two variable 
case (double zeta values) was already studied by Euler. 

MZVs in the general case were introduced 1990s with motivations from number theory IT701 , 
combinatorics [|40l and quantum field theory IfTTTl . Since then the subject has turned into an active 
area of research that involves many areas of mathematics and mathematical physics Ifi31l . Its 
number theoretic significance can be seen from the fact that all MZVs are periods of mixed Tate 
motives over Z and the conjecture that all periods of mixed Tate motives are rational combinations 
of MZVs [|25l|23[69l. 

It has been discovered that the analytically defined MZVs satisfy many algebraic relations. Fur- 
ther it is conjectured that these algebraic relations all follow from the combination of two algebra 
structures: the shuffle relation and the stuffie (harmonic shuffle or quasi- shuffle) relation ||46l . 
This remarkable conjecture not only links the analytic study of MZVs to the algebraic study of 
double shuffle relations, but also implies the more well-known conjecture on the algebraic inde- 
pendence of £(2), £(2k +l),k>\, over Q. 
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Many results on algebraic relations among MZVs can be regarded as generalizations of Eu- 
ler's sum formula and decomposition formula on double zeta values which preceded the general 
developments of multiple zeta values by over two hundred years. We summarize these results 
in Section [3] With the non-experts in mind, we first give in Section |2] preliminary concepts and 
results on double shuffle relations for MZVs and the related Rota-Baxter algebras. 

1.2. Renormalization. Values of the Riemann zeta function at negative integers are defined by 
analytic continuation and possess significant number theory properties (Bernoulli numbers, Rum- 
mer congruences, p-adic L-functions, • • • ). Thus it would be interesting to similarly study MZVs 
outside of the convergent domain of the corresponding nested sums. However, most of the MZVs 
remain undefined even after the analytic continuation. To bring new ideas into the study, we 
introduce the method of renormalization from quantum field theory. 

Renormalization is a process motivated by physical insight to extract finite values from diver- 
gent Feynman integrals in quantum field theory, after adding in a so-called counter-term. Despite 
its great success in physics, this process was well-known for its lack of a solid mathematical 
foundation until the seminal work of Connes and Rreimer lfT4l[T5l[T6l[50ll . They obtained a Hopf 
algebra structure on Feynman graphs and showed that the separation of Feynman integrals into 
the renormalized values and the counter-terms comes from their algebraic Birkhoff decomposition 
similar to the Birkhoff decomposition of a loop map. 

The work of Connes and Rreimer establishes a bridge that allows an exchange of ideas between 
physics and mathematics. In one direction, their work provides the renormalization of quantum 
field theory with a mathematical foundation which was previously missing, opening the door 
to further mathematical understanding of renormalization. For example, the related Riemann- 
Hilbert correspondence and motivic Galois groups were studied by Connes and Marcolli IfTTll . 
and motivic properties of Feynman graphs and integrals were studied by Bloch, Esnault and 
Rreimer [5J. See [|2l[IU[56l for more recent studies on the motivic aspect of Feynman rules and 
renormalization. 

In the other direction, the mathematical formulation of renormalization provided by the alge- 
braic Birkhoff decomposition allows the method of renormalization dealing with divergent Feyn- 
man integrals in physics to be applied to divergent problems in mathematics that could not be dealt 
with in the past, such as the divergence in multiple zeta values ||36l 1371 1731 1551 and Chen symbol 
integrals Il54ll55ll . We survey these studies on renormalization in mathematics in Sections [5J and 
[6] after reviewing in Section |4] the general framework of algebraic Birkhoff decomposition in the 
context of Rota-Baxter algebras. We further present an alternative renormalization method using 
Speer's generalized evaluators ll67l and show it leads to the same renormalized double zeta values 
as the algebraic Birkhoff decomposition method. 

We hope our paper will expose this active area to a wide range of audience and promote its 
further study, to gain a more thorough understanding of the double shuffle relations for convergent 
MZVs and to establish a systematical renormalization theory for the divergent MZVs. One topic 
that we find of interest is to compare the various renormalization methods presented in this paper 
from an abstract point of view in terms of a renormalization group yet to be described in this 
context, again motivated by the study in quantum field theory. With implications back to physics 
in mind, we note that MZVs offer a relatively handy and tractable field of experiment for such 
issues when compared with the very complicated Feynman integral computations. 
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2. Double shuffle relations for convergent multiple zeta values 

All rings and algebras in this paper are assumed to be unitary unless otherwise specified. Let k 
be a commutative ring whose identity is denoted by 1 . 

2.1. Rota-Baxter algebras. Let A e k be fixed. A unitary (resp. nonunitary) Rota-Baxter k- 
algebra (RBA) of weight A is a pair (R, P) in which R is a unitary (resp. nonunitary) k-algebra 
and P : R — » R is a k-linear map such that 

(1) P(x)P(y) = P(xP(y)) + P(P(x)y) + AP(xy), Vx, y e R. 

In some references such as Il55ll , the notation 6 = -A is used. 

We will mainly consider the following Rota-Baxter operators in this paper. See lfP9ll30l[64ll for 
other examples. 

Example 2.1. (The integration operator) Define the integration operator 



(2) /(/)(*)= f(t)dt 



-f 

Jo 



on the algebra C[0, oo) of continuous functions f(x) on [0, oo). Then it follows from the integra- 
tion by parts formula that / is a Rota-Baxter operator of weight @|. 

Example 2.2. (The summation operator) Consider the summation operator 11751 

P(f)(x):=J]f(x + n). 

Under certain convergency conditions, such as f(x) = 0(x~ 2 ) and g(x) = 0(x~ 2 ), P(f)(x) and 
P(g)(x) define absolutely convergent series and we have 

P(f)(x)P(g)(x) = Y J f(x + m)Y J 8(x + n) 

(3) = ^ f(x + m)g(x + n)+ ^ f(x + m)g(x + n) + ^ f(x + m)g(x + m) 

= P(fP(g))(x) + P(gP(f))(x) + P(fg)(x). 
Thus P is a Rota-Baxter operator of weight 1. 

Example 2.3. (The partial sum operator) The operator P defined on sequences cr : N — > C by: 

n 

(4) P(cr)(n) = Y J o-(k) 

k=0 

satisfies the Rota-Baxter relation with weight -1. Similarly, the operator Q = P - Id which acts 
on sequences cr : N — » C by: 

n-\ 

(5) Q(cr)(n) = Y J cr(k) 

k=Q 
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satisfies the Rota-Baxter relation with weight 1 . 

Example 2.4. (Laurent series) Let A = k[s~ l ,s\\ be the algebra of Laurent series. Define 
n : A -> A by 

nC^OaO = J]a n s n . 

n n<0 

Then IT is a Rota-Baxter operator of weight -1. 

2.2. Shuffles, quasi-shuffles and mixable shuffles. We briefly recall the construction of shuffle, 
stuffle and quasi- shuffle products in the framework of mixable shuffle algebras ll32ll33Tl . 

Let k be a commutative ring. Let A be a commutative k-algebra that is not necessarily unitary. 
For a given A e k, the mixable shuffle algebra of weight A generated by A (with coefficients in 
k) is MS (A) = MSk, / i(A) whose underlying k-module is that of the tensor algebra 

(6) T(A) = (J)A®* = keA©A® 2 ©--- 

k>0 

equipped with the mixable shuffle product o /{ of weight A defined as follows. 

For pure tensors a = a\ ® . . . <8> a m 6 A® m and b = b\ ® . . . <g> b n e A®", a shuffle of a and b is 
a tensor list of a, and bj without change the natural orders of the a,-s and the bjS. More precisely, 
for o~ G £ w := {r e 5^ | t _1 (1) < ••• < t" 1 ^), t" 1 ^ + 1) < ••• < T~\k + f)}, the shuffle of a and 
b by cr is 

o ffl[r b := Ccr(1) • • • 9 Co-(k+t) > where Cj = | ^ ^ { <^ * + ^ 
The shuffle product of a and b is 

dob := ^ auio-b. 

More generally, for a fixed A e k, a mixable shuffle (of weight /I) of a and b is a shuffle of 
a and b in which some (or none) of the pairs a, <S> bj are merged into A a\bj. Then the mixable 
shuffle product of weight A is defined by 

(7) ao ,ib = ^ mixable shuffles of a and b 

where the subscript A is often suppressed when there is no danger of confusion. For example, 

aiO^Z?! (8> b-i) := a\ ® b\ <8> b-i + b\ ® a\ ® bi + b\ ® Z?2 ® «i + ^i(ai^i) ® Z7 2 + ^£>i <8» {a\b-i) . 

" ^ ' v ' 

shuffles merged shuffles 

With 1 e k as the unit, this product makes T(A) into a commutative k-algebra that we denote by 
MS k /! (A). See [|32l for further details on the mixable shuffle product. When A = 0, we simply 
have the shuffle product which is also defined when A is only a k-module, treated as an algebra 
with the zero multiplication. 

We have the following relation between mixable shuffle product and free commutative Rota- 
Baxter algebras. A Rota-Baxter algebra homomorphism / : (R, P) — » (R',P') between Rota- 
Baxter k-algebras (R, P) and (R', P') is a k-algebra homomorphism / : R — » R' such that / o P = 
P'of. 
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Theorem 2.5. The tensor product algebra 111(A) := HI M (A) = A <g> MS kj/i (A), wzY/z the 

linear operator P A : 111(A) — > 111(A) sending a — > 1 ® a, w /ree commutative Rota-Baxter 
algebra generated by A in the following sense. Let j& : A —* 111(A) be the canonical inclusion 
map. Then for any Rota-Baxter k-algebra (R, P) and any k-algebra homomorphism <p : A — > R, 
there exists a unique Rota-Baxter k-algebra homomorphism <p : (111(A), P A ) — » (R, P) such that 
<p = <p o j A as k-algebra homomorphisms. 

The product o A can also be defined by the following recursion lfT8l [371 15511 which provides the 
connection between mixable shuffle algebras and quasi- shuffle algebras of Hoffman [|42l . First 
we define the multiplication by A 80 = k to be the scalar product. In particular, 1 is the identity. 
For any m, n > 1 and a := a\ <S> ■ ■ ■ <8> a m e A®'", b := &i ® • • • <S> b„ e A®", define ao\b by induction 
on the sum m + n > 2. When m + = 2, we have a = a\ and b = b\. We define 



Assume that cto^b has been defined for m + n > > 2 and consider a and b with m + n = k + 1 . 
Then m + /O 3 and so at least one of m and n is greater than 1 . We define 



' ai®bi®---®b„ + bi® (aiO A (b 2 ® • • • ® &„)) 

+/l(aiZ?i) ® b 2 <S> • • • ®b n , when m = l,n > 2, 
«i ® ((«2 ® ■ • ■ ® a m )o,\bi) + b\ <8> «i <g> • • • ® a m 

+/l(a[Z7i) ® a 2 <S> • • • <S> a m , when m > 2, n = 1, 
ai <g> ((a 2 <8> • • • ® a m )o A {bi ® • • • <g> />„)) + £i <g> ((ai <g> • • • <S> a m )o A (b 2 ® • • • ® b n )) 

+A{a\b\){{a 2 <S> ■ • ■ <8> a m )o A (b 2 <S> • • • <8> £„)), when m, zO 2. 



Here the products by o /{ on the right hand side of the equation are well-defined by the induction 
hypothesis. 

Let S be a semigroup and let kS = ^Zses k 5 De me semigroup nonunitary k-algebra. A canon- 
ical k-basis of (kSf k , k>0,is the set S® k := {ji <g> • • • ® s k \ s t e S, 1 < i < /<:}. Let 5 be a graded 
semigroup 5 = JJ i>0 S„ S,5 7 c S i+ j such that < oo, z > 0. Then the mixable shuffle product 
<>! of weight 1 is identified with the quasi-shuffle product * defined by Hoffman 11421 [T8l |3~7i 

Notation 2.6. To simplify the notation and to be consistent with the usual notations in the litera- 
ture on multiple zeta values, we will identify Si <8> • • • <S> Sk with the concatenation s\ • • ■ Sk unless 
there is a danger of confusion. We also denote the weight 1 mixable shuffle product o { by * and 
denote the corresponding mixable algebra MSk,i(A) by !K* A . Similarly, when A is taken to be a 
k-module, we denote the weight zero mixable shuffle algebra MS ki0 (A) by "K A . 

Yet another interpretation of the mixable shuffle or quasi-shuffle product can be given in terms 
of order preserving maps that are called stuffle in the study of MZVs but could be traced back to 
Carrier's work [[T2l on free commutative Rota-Baxter algebras. 

For positive integers k and I, denote [k] = {1, • • • ,k] and [k + 1 , k + €\ = [k + 1, • • • ,k + £}. 
Define 



Let a e A® , b e A® and {up, if/) e We define am^ ^b to be the tensor whose z'-th factor is 



ao A b = a\ ® b\ + b\ <8> a\ + \a\b\. 



(8) 




cp : [k] — » [k + €], if/ : [■£] — > [k + £] are order preserving 
injective maps and im(^) U im(i/0 = [k + €\ 



(9) 



( 



aj if i = ip(j) 
bj if i = i/f(J) 



1 < i < k + 1, 
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with the convention that a% = b% = 1 . Then we have 

(10) ctmb = ^ ctm^b. 



More generally, for < r < min(/t, €), define 



cp : [k] — > [k + (, - r], ft : [€] — » [k + 1 - r] are order preserving 
injective maps and im(^) U im(i/r) = [k + t - r] 

Clearly, J^/.o = ^u- Let o e A® k , b e A® f and (<p, if/) e "Jk,e, r - We define Cta^b to be the tensor 
whose i-th factor is 



dj if i = (p(j), i i'mwjj 
bj if i = if/(j), i£im<p \ = VHoV 1 ©' 1 < z ' < k + ^ _ r ' 



a^v if z = <p(j), i = if,(f) 
with the convention that a = = 1 . Then we have [|32] [361 

min(it,f) 

(11) ao /t b= ^ i r ( ^ am^b). 
In particular, 

o * b = ^ ( Yj am (v^) b ) = X am ^ b 

where V = 

Equivalently, let stfl(/t, r) denote the set of surjective maps from [k + £] to [k + £ - r] that 
preserve the natural orders of [k] and {k + 1, • • • + Let 

stfl(/c, = [J stfl(/cj,r). 

r=0 

Then 

(12) (fli ® • • • ® a*) * (a fc+1 ® • • • ® a* +f ) = ^ cf ® • • • ® c^_ r , cf = a,-. 

A connected filtered Hopf algebra is a Hopf algebra (H, A) with k-submodules H (n) , n > 
of H such that 

#(«) c H in+l) , U n>Q H (n) = H, H^H® c H ( p +c '\ 
A(H {n) ) c 2 p+ 4=„ fl0») ® #C«), H (0) = k (connectedness). 

On the algebra MS k /! (A) further define 

A(ai ® • • • ® a n ) = 1 (S?)(ai ® • • • ®a„) + a x <S^)(a 2 ® • • • ® a n ) 



(13) +••• + (ai ® ■ ■■a n -i)QQa n + (a\ ® • • • ®a„) ® 1. 

Then A extends by linearity to a linear map MS kj/i (A) -> MS M (A) MS ki/( (A). 
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Theorem 2.7. ([[37l|42l|55lJ The triple (MS ki/l (A), o A , A), together with the unitu : k ^ MS M (A) 
and the counit s : MS ki/ i(A) — > k projecting onto the direct summand k c MS ki/ i(A), equips 
MSk ^(A) wzY/z structure of a connected filtered Hopf algebra with the filtration MS(A) ( " ) := 
V 4®' 

We also have the following easy extension of Hoffman's isomorphism between the shuffle Hopf 
algebra and the quasi-shuffle Hopf algebra (see also |[T8T0 . Recall the notation "K A = MS ki i(A) 
and IK™ = MS k , (A). 

Theorem 2.8. ([|42l[55l) Let kbea Q-algebra. There is an isomorphism of Hopf algebras : 

(14) e\p:^C^JC A . 

Hoffman's isomorphism (fl4)) is built explicitly as follows. Let CP(n) be the set of compositions 

of the integer n, i.e. the set of sequences / = (z'i , . . . , 4) of positive integers such that i\ -\ h4 = n. 

For any u = v\ ® • • • ® v„ e T(A) and any composition / = . . . , it) of n we set: 

I[u] := (vi v h ) ® (v, 1+1 v, 1+i2 ) ® • • • ® (v il+ ... +4 _ 1+1 v„). 

Then the isomorphism exp is defined by 

expw = V . , 1 . I[u]. 

^— ' ^ *1 ' ' ' ' Ik- 

Moreover (021, Lemma 2.4), the inverse log of exp is given by : 

\ogu= ) -/[«]. 

2.3. Double shuffle of MZVs and related conjectures. A multiple zeta value (MZV) is defined 
to be 



y — — 



(15) £(su--- ,s k ) :■■ / . Sk 

t—L , n, ■ ■ -n, 

where s, > 1 and s\ > 1 are integers. As is well-known, an MZV has an integral representation 
due to Kontsevich ||5T| 



J J 

Jo Jo 



/i(*i) f\4$\A) 
J 1 - tj, j = S U Si + S 2 , ■ ■ ■ , 51 + • • • + S k , 



(i6) a*u---,sk) 

Here |.s] = si + • • • + Sk and 

W> - l ,„ otherwise. 

The MZVs spanned the following Q-subspace of R 

MZV := Qiasi, ■ ■ ■ , s k ) I s t > 1, s l > 2} c R. 

Since the summation operator in Eq. (TT3T ) and the integral operator in Eq. © are both Rota- 
Baxter operators (of weight 1 and respectively) by Example 12.21 and Example 12. 1[ it can be 
expected that the multiplication of two MZVs follows the multiplication rule in a free Rota-Baxter 
algebra and thus in a mixable shuffle algebra. This viewpoint naturally leads to the following 
double shuffle relations of MZVs. 
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For the sum representation of MZVs in Eq. (TT5T) . consider the semigroup 

Z := {z s | s e Z>i, z s • z t = z s+t , s,t> 1.} 

With the convention in Notation 12.61 we denote the quasi-shuffie algebra "K* := IKqz which 
contains the subalgebra 

JHS:=Q©(0Qz„ ••■**)■ 

Sl>\ 

Then the multiplication rule of MZVs according to their summation representation follows from 
the fact that the linear map 

(17) C ■ ^ -» MZV, z,,,.,,, h-> • • • , s k ) 

is an algebra homomorphism []4Tll46l . 

For the integral representation of MZVs in Eq. (fT6l ), consider the set X = {x , X\). With the 
convention in Notation l2.6[ we denote the shuffle algebra IK m := IKqx which contains subalgebras 

Then the multiplication rule of MZVs according to their integral representations follows from the 
statement that the linear map 

r : IK m ^ MZV, xj 1 - 1 ^ • • • < _1 xi i-> {(su ■ • ■ , J*) 

is an algebra homomorphism [|4Tll46l . 

There is a natural bijection of Q- vector spaces (but not algebras) 

77 : K^ IK*, 1 <-» 1, • ■■Jtg t_1 j:i «-> z„ t ._,, t . 

that restricts to a bijection of vector spaces rj : IKo 1 — » IKo. Then the fact that MZVs can be 
multiplied in two ways is reflected by the commutative diagram 

<rsm 1 ^ or* 




MZV 

Through 77, the shuffle product m on IKf and IKo 1 transports to a product m, on IK* and IKq. That 
is, for wi,w 2 6 IKq, define 

(18) wimw 2 := r](rf l {w\)mrf l (w 2 )). 
Then the double shuffle relation is simply the set 

{wim, w 2 - W\ * w 2 I W\ , w 2 s IKo*} 
and the extended double shuffle relation (461 [63l |75]| is the set 

(19) {wim w 2 - W\ * w 2 , zim w 2 - Zi * w 2 1 Wi, w 2 e IKo}. 

Theorem 2.9. ( ll4Tll46ll63lO Le? Teds ^ ^ idea/ o/IKq generated by the extended double shuffle 
relation in Eq. ( fTPl) . 77ien Teds z * zm ^ kernel of '£* '. 

It is conjectured that 7 E ds is in fact the kernel of £*. A consequence of this conjecture is the 
irrationality of £(2n + 1), n > 1. 
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3. Generalizations of Euler's formulas 

We begin with stating Euler's sum and decomposition formulas in Section l3TTl Generalizations 
of Euler's sum formula are presented in Section [321 and generalizations of Euler's decomposition 
formula are presented in Section! 



3.1. Euler's sum and decomposition formulas. Over two hundred years before the general 
study of multiple zeta values was started in the 1990s, Goldback and Euler had already considered 
the two variable case, the double zeta values ll23l l66ll 

£(sus 2 ):= V. — 3r-ia ■ 

n\>ii2>l 1 1 

Among Euler's major discoveries on double zeta values are his sum formula 

n-\ 

{(i, n-i) = £(n) 



i=2 



expressing one Riemann zeta values as a sum of double zeta values and the decomposition for- 
mula 

(20) S(r)as) = ^ W + *, j - *) + £ I"*' 1 )&s + k,r- k), r,s>2, 

k=0 ' ' k=0 * ' 

expressing the product of two Riemann zeta values as a sum of double zeta values. 

A major aspect of the study of MZVs is to find algebraic and linear relations among MZVs, 
such as Euler's formulas. Indeed a large part of this study can be viewed as generalizations of 
Euler's formulas. 

3.2. Generalizations of Euler's sum formula. Soon after MZVs were introduced, Euler's sum 
formula was generalized to MZVs (140112811711 as the well-known sum formula, followed by quite 
a few other generalizations that we will next summarize. 

3.2.1. Sum formula. The first generalization of Euler's sum formula is the sum formula conjec- 
tured in [gOl. Let 

(21) I(n,k) = {(*!, ■•• ,s k ) | Si + ■ ■ ■ + s k = n,Si > l,s 1 > 2}. 

For ,?= (si, ■ ■ ■ , Sk) 6 I(n, k), define the multiple zeta star value (or non-strict MZV) 

(22) X ,, 1 r k - 

n\>~>rik>\ 1 k 

Note the subtle different between the notations in Eq. (fTTb and [* in Eq. (|22j) . 
Theorem 3.1. (Sum formula) For positive integers k < nwe have 

(23) £ = E ^ = (l-l)^ (n) - 

sel(n,k) kel(n,k) * ' 
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The case of k = 3 was proved by M. Hoffman and C. Moen 11431 and the general case was proved 
by Zagier 11711 with another proof given by Granville ll2~8Tl . Later S. Kanemitsu, Y. Tanigawa, M. 
Yoshimoto ||48T| gave a proof for the case of k = 2 using Mellin transformation. 

J.-I. Okuda and K. Ueno ll62ll gave the following version of the sum formula 



£(::)(£«*) =(:>»> 

l~r \ ' fcTIn 1A \ I 



k=r x ' sel(n,k) 

for n > r > 1 from which they deduced the sum formula Eq. (l23l) . 

3.2.2. Ohno's generalized duality theorem. Another formula conjectured in [|40ll is the duality 
formula. To state the duality formula, we need an involution r on the set of finite sequences of 
positive integers whose first element is greater than 1 . If 

?=(1+& 1 ,1 1 —_ 1 1,-- - A + b^h—^l), 

then 

r(s) = (1 + a k , 1,- - • ,1, ■••,!+ ai,l," • • , !)• 

b k -\ fei-i 

Theorem 3.2. (Duality formula) 

= f(r(^)). 



This formula is an immediate consequence of the integral representation in Eq. (fT6l) . 
Y. Ohno 11571 provided a generalization of both the sum formula and the duality formula. 



Theorem 3.3. (Generalized Duality Formula (57J) For any index set s = (si, ■ ■ ■ , s k ) with 
Si > 2, S2 > 1, • • • , Sk > 1, and a nonnegative integer I, set 

Z(s u ■ ■ ■ , s k ;€) = ^ g(si + c u -'- ,s k + c k ). 

C\ + ■ • ■ + Q = t 

Ci>0 

Then 

Z(s*;€) = Z(r(s)-J). 

When t = 0, this is just the duality formula. When s = (k + 1) and I = n — k — 1, this becomes 
the sum formula. 

3.2.3. Sum formulas with further conditions on the variables. M. Hoffman and Y. Ohno p4| 
gave a cyclic generalization of the sum formula. 

Theorem 3.4. (Cyclic sum formula) For any positive integers s\, ■ ■ ■ ,s k with some Sj > 2, 

k s-2 

2^{(Si+ l,s i+ i,'-- ,s k ,s u --- ,Jf-i) = z\j ~ j' Si+U ' " ' s k> s w- >Si- U j+ !)• 

i=l {i\si>2} j=0 

Y. Ohno and N. Wakabayashi ll59ll gave a cyclic sum formula for non-strict MZVs and used it 
to prove the sum formula Eq. (T23l) . 
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Theorem 3.5. (Cyclic sum formula in the non-strict case) For positive integers k < n and 
(s\, ■ ■ ■ , s k ) G I(n, k) we have 



k s-2 



(24) 



Y Y £*(Si - j, s M , ■ ■ ■ , s k , s u ■ ■ ■ , Si-u j + 1) = n£(n + 1), 

i=l j=0 



where the empty sums are zero. 

M. Eie, W.-C. Liaw and Y. L. Ong ||22l| gave a generalization of the sum formula by allowing 
a more general form in the arguments in the MZVs. 

Theorem 3.6. For all positive integers n, k with n > k, and a nonnegative integer p, 



E 



£(c\, ■ ■ ■ ,c p+ i) 



s\ + ■ ■ ■ + s/c = n 

si>2 



ci+ ■■■ + Cp+i =n+p 
c\ > n — k+l 



When p = 0, this becomes the sum formula. 

Y. Ohno and D. Zagier ll60l studied another kind of sum with certain restrictive conditions. Let 
I(n,k,r) = {(s u ■■■ ,s k )\si£ Z >u s x + ■ ■ ■ + s k = n,#{s t \ Si > 2} = r} 

and put 

G(n,k,r)= fCfr 

fel(n,k,r) 

They studied the associated generating function 

o>(x, y ,z)= G{n - K r )* n ~ k ~ r y k ~ r z r ~ l e y. z] 



r>l,k>r,n>k+r 



and proved the following 
Theorem 3.7. We have 



1 



/ CO 



1 - exp 



xy - z 

where S n (x,y,z) are given by the identity 

xy - z 



\n=2 



> . S n (x,y,z) 

n 



log 1 



(l-x)(l-y) 



Z 

n=2 



S n (x,y,z) 



and the requirement that S n (x,y,z 2 ) is a homogeneous polynomial of degree n. In particular, 
all of the coefficients G(n,k,r) can be expressed as polynomials in £(2), f(3), ... with rational 
coefficients. 
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3.2.4. Sum formulas for q-MZVs. The concept of ^-multiple zeta values (g-MZVs, or multiple 
g-zeta values) was introduced as a "quantumization" of MZVs that recovers MZVs when q i-» 
1 [01221. 

For positive integers S\, • • ■ ,s k with s\ > 2, define the q-MZV 

„ni(il-l)+-+n*(st-l) 



and the non- strict q-MZV 



„ni(il-l)+-+n t (it-l) 



[m]*i • • ■ 

where [«] = -t- 2 -. 

D. M. Bradley [9] proved the ^-analogue of the sum formula for £ q . 

Theorem 3.8. (^-analogue of the sum formula) For positive integers < k < nwe have 

(25) 2_j £«( f i> " ' " > s k) = £?(")■ 

S] H hi|t=« 

Y. Ohno and J. -I. Okuda [58] gave the following ^-analogue of the cyclic sum formula (1241) and 
then used it to prove a ^-analogue of the sum formula for £*. 

Theorem 3.9. (^-analogue of the cyclic sum formula) For positive integers < k < n and 
(si, ■ ■ ■ , Sk) G I(n, k) we have 

k Sj—2 k j \ 

i=l 7=0 f=0 > ' 

where the empty sums are zero. 

Theorem 3.10. (^-analogue of the sum formula in the non-strict case) For positive integers 
< k < n we have 

1 / - \ k ~ l 

(26) ■•■•'*> = — rti E (n_1 - " " ^ 

Sl+-+iJt=« 

3.2.5. Weighted sum formulas. In the other direction to generalize Euler's sum formula, there is 
the weighted version of Euler's sum formula recently obtained by Ohno and Zudilin |f6TTl . 

Theorem 3.11. (Weighted Euler's sum formula flUI) For any integer n > 2, we have 

n-l 

(27) Y i 2 i ai,n-i) = (n + l)£(ri). 

i=2 

They applied it to study multiple zeta star values. By the sum formula, Eq. (1271) is equivalent 
to the following equation 

n-l 

(28) Yj(2 i -l)ai,n-i) = n^n). 

i=2 
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As a generalization of Eq. (1281) . two of the authors proved the following 
Theorem 3.12. (Weighted sum formula 11351 ) For positive integers k>2 and n > k+l, we have 



k-\ 



£ + - 1)(( ^ 2 St -°^ 1) ) + 2 s *— (k - z) )]{(su ■••,**) = n{(n), 

Sj>l,si>2 i'=2 
S\ H h^=n 

where 5; = si + h j = 1, • • • , & - 1. 

3.3. Generalizations of Euler's decomposition formula. Unlike the numerous generalizations 
of Euler's sum formula, no generalization of Euler's decomposition formula to MZVs, neither 
proved nor conjectured, had been given until |T34ll even though Euler's decomposition formula 
was recently revisited in connection with modular forms Il24ll and weighted sum formula [1611 on 
weighted sum formula of double zeta values, and was generalized to the product of two g-zeta 
values [[T01I72I. 

3.3.1. Euler's decomposition formula and double shuffle. The first step in generalizing Euler's 
decomposition formula is to place it as a special case in a suitable broader context. In [341 Euler's 
decomposition formula was shown to be a special case of the double shuffle relation. We give a 
proof of Euler's formula in this context before presenting its generalization in the next subsection. 
We recall that the extended double shuffle relation is the set 

{wim w 2 - W\ * w 2 , Zim w 2 -zi*w 2 \ W\,w 2 e %o). 

Thus the determination of the double shuffle relation amounts to computing the two products * 
and m . 

It is straightforward to compute the product *, either from its recursive definition in Eq. © 
or its explicit interpretation as mixable shuffles in Eq. © and stuffles in Eq. (fTTT) or (PT21) . For 
example, to determine the double shuffle relation from multiplying two Riemann zeta values £(r) 
and £(s), r, s > 2, one uses their sum representations and easily gets the quasi-shuffle relation 

On the other hand, computing the product m is more involved as can already be seen from 
its definition in Eq. (fl"8l) . One first needs to use their integral representations to express £(r) and 
£(s) as iterated integrals of dimensions r and s, respectively. One then uses the shuffle relation to 

express the product of these two iterated integrals as a sum of y j iterated integrals of dimension 

r + s. Finally, these last iterated integrals are translated back to MZVs and give the shuffle relation 
of £(r)£(s). As an illustrating example, consider £(100)^(200). The quasi-shuffle relation is 
simply £(100)^(200) = £(100, 200) + £(200, 100) + £(300), but the shuffle relation is a large sum 

of ^ ioo ) snu ^ es °f l en S m (dimension) 300. As we will show below, an explicit formula for this 

is precisely Euler's decomposition formula (1201) . 

Theorem 3.13. For r,s>2, we have 

r-1 



Zt r+k-l\ ( s+k-l\ 

I I Zr+kZ s -k + 2_j\ k I Z s+kZr-k- 

k=Q \ I k=Q ' ' 



14 LI GUO, SYLVIE PAYCHA, BINGYONG XIE, AND BIN ZHANG 

Via the algebra homomorphism in Eq. (fTTl) this theorem immediately gives Euler's decom- 
position formula. Applying to the above example, we have 

199 / v 99 i v 

£(100)£(200) = Yj 100 ^ 1 £( 100 + k > 200 - k "> + Y\ 2 °T 1 ^ 20 ° + k ' 100 ~ k) - 

k=o ' I k=0 ' ' 

Proof. Following the definition of m, in Eq. (PT8T) . we have 

So we just need to prove 

Xq 1 XiuiXq l X\ = t )"' C +A 1 - V:iX U + ^ I . \ X +k ^l^O * ^l 

k=0 ' ' <fc=0 > ' 

since m (x^^xiXg - * -1 *!) = z r +kZ s -k and m (x^* -1 *^ - * -1 .*!) = z s +kZ r -k- This has a direct shuffle 
proof JU. But we use the description of order preserving maps of shuffles in order to motivate the 
general case. 

By Eq. (fTOl) . we have 

r- 1 s-l r-1 s-l 

X XiiiiXq Xi — ?, x x l m {if,ifi) x Q x 0- 

Since tp and if/ are order preserving, we have the disjoint union 3(r, s) = J(r, s)' U J(r, 5)" where 

J(r,J)' = {(p^)e3(r,i)|^(s) = r + f] 

and 

3(r,s)" = {(cp,if/)e3(r,s)\(p(r) = r + s}. 

Again by the order preserving property, for ((p, if/) e J(r, 5)', we must have <p(r) = r + k where 
k > 0. Thus for such (<£, t/0, we have 

Xq Xiui^^Xq X\ = Xq X\Xq X\ 

since im ip U im ^ = [r + 5] . For fixed > 0, ^>(r) = r + k means that there are k elements i\, ■ ■ • ,4 
from [s - 1] such that if/(ij) £ [r + k — 1] since ^(s) = r + s. Thus k > 5 - 1 and, since t/r is order 

(r+k-l \ 
^ I such i/f's since <^([fc]) can take any 

k places in [r + k - 1] in increasing order and then (f>([r]) takes the rest places in increase order. 
Thus 

Zr-l s-l X" 1 ( ' r+k-l\ r+k-\ s-Jfc-1 

•^0 ^l™^)-^ X l = I , \ X X 1 X X l- 
(if,t/i)eJ(r,s)' k=0 * ' 

By a similar argument, we have 

r-l 



^ *o lx i m (^^) x o ljc i - Yj\ k ) 



.v+A-1 r-£-l 

Aq A J -A. Q -A. J • 



(p,^)63(r,i)" k=Q 

This completes the proof. 
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3.3.2. Generalizations of Euler's decomposition formula. In a recent work [34], two of the au- 
thors generalized Euler's decomposition formula in two directions, from the product of one vari- 
able functions to that of multiple variables and from multiple zeta values to multiple polyloga- 
rithms. 

A multiple polylogarithm value [|71|25l|26l is defined by 



where \z,-\ < 1, s t e Z >{ , 1 < i < k, and (*i,Zi) + (1, 1). When zi = 1,1 < i < k, we obtain the 
multiple zeta values £(si, • • • , s k ). With the notation of 0, we have 



where (bu--- ,b k ) = (z^AziZi) \ ■ ■ ■ , (z\ • • -Zk) l ). 

To state the result, let k and t be positive integers and let J k j be as defined in Eq. ©. Let 
t = (r u ■ ■ ■ , r k ) € Z k >v f= (s u ■ ■ ■ , s ( ) 6 and t = (t u ■ ■ ■ , t kU ) e Z* + / with \r) + \s] = Here 

\r\ = r x -\ + r k and similarly for |^ and \t\. Denote R { = r\ H + r t for i e [£], 5, = Si H + j,- 

for z 6 [£] and T { = t\ + ■ ■ ■ + t t for i e [k + £]. For ((p, iff) 6 J w and ie[Hf], define 




Liji,-,*t(zi. ' ' ' ,Zk) = K 



s\,--- ,s k 
bu--- ,h 



>=- I 



V^/ \b 2 ' • ■ ■ V b k > 

S I So Sl- 

n l n 2'-- n k 



i = h 




with the convention that r% = s® = 1 . 
With these notations, we define 



(29) 




For £ e (S 1 ) t and £ e (5 as in Eq. ©, define 

(30) am (lf ^)b = («^-l (1)^-1 (1), • • • , «^-l(fc+f)^-l(A:+f))- 

Theorem 3.14. ( Il34l ) Lef /c, ^ £e positive integers. Let r e Z* j and se Z> r Le? a = (ai, • ■ • ,a k ) € 
(S l ) k andt = (b u ■ ■ ■ , b e ) e (S l ) c such that [ n ] ± [ 1 ] and [ 11 ] * [ 1 ]. 77ien 





where c^\i) is given in Eq. ([29]) and dm^^b is given in Eq. 
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Corollary 3.15. Let r e Z*j and s€ Z^j with r\, si > 2. Then 

k+e 

r*€Z*^,|rf=|?l+|j1 (v>,<l>)£\t i=l 

where c^'^ii) is given in Eq. ri29l) . 

4. The algebraic framework of Connes and Kreimer on renormalization 

The Algebraic Birkhoff Decomposition of Connes and Kreimer is a fundamental result in their 
ground breaking work |[T5l on Hopf algebra approach to renormalization of perturbative quantum 
field theory (pQFT). This decomposition also links the physics theory of renormalization to Rota- 
Baxter algebra that has evolved in parallel to the development of QFT renormalization for several 
decades. 

The introduction of Rota-Baxter algebra by G. Baxter 01 in 1960 was motivated by Spitzer's 
identity [68] that appeared in 1956 and was regarded as a remarkable formula in the fluctuation 
theory of probability. Soon Atkinson [3] proved a simple yet useful factorization theorem in 
Rota-Baxter algebras. The identity of Spitzer took its algebraic form through the work of Cartier, 
Rota and Smith 021 [65l (1972). 

It was during the same period when the renormalization theory of pQFT was developed, 
through the the work of Bogoliubov and Parasiuk [6] (1957), Hepp [[391 (1966) and Zimmer- 
mann [1741 (1969), later known as the BPHZ prescription. 

Recently QFT renormalization and Rota-Baxter algebra are tied together through the algebraic 
formulation of Connes and Kreimer for the former and a generalization of classical results on 
Rota-Baxter algebras in the latter [|20l |2~H . More precisely, generalizations of Spitzer's identity 
and Atkinson factorization give the twisted antipode formula and the algebraic Birkhoff decom- 
position in the work of Connes and Kreimer. 

We recall the algebraic Birkhoff decomposition in Section |4~T] prove the Atkinson factorization 
in Section 14.21 and derive the algebraic Birkhoff decomposition from the Atkinson factorization 
in Section 1431 

4.1. Algebraic Birkhoff decomposition. For a k-algebra A and a k-coalgebra C, we define the 
convolution of two linear maps /, g in Hom(C, A) to be the map / ★ g 6 Hom(C, A) given by the 
composition 

A f®g m 

C — > C <S> C — > A <g> A — > A. 

Theorem 4.1. (Algebraic Birkhoff Decomposition) Let H be a connected filtered Hopf algebra 
over C. Let (A, II) be a commutative Rota-Baxter algebra of weight -1 with n 2 = II. 

(a) For (f> e char(H, A), there are unique linear maps : H — > k + 11(A) and <p + : H — » 
k + (id - II)(A) such that 

(31) = 0! ( - 1) *0 + . 
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(b) The elements cf)- and </>+ take the following forms on kere. 

(32) 0_(x) = -U(<f>(x) + £ 0_(x')0(*")), 

(33) + (x) = ft(0(x) + J] 0-(*W')), 

where we have used the notation A(x) = 1 <g> jc + <8> 1 + 2 m ® with x', x" e ker s. 

(c) 77ze linear maps </>- anJ 0+ are a/50 algebra homomorphisms. 

We call + the renormalization of (f> and call </>_ the counter-term. Here is roughly how the 
renormalization method can be applied through the Algebraic Birkhoff Decomposition. See the 
tutorial article OTI for further details, examples and references. 

Theorem 14.11 can applied to renormalization as follows. Suppose there is a set of divergent 
formal expressions, such as MZVs with not necessarily positive arguments, that carries a certain 
algebraic combinatorial structure and from which we would like to extract finite values. On one 
hand, we first apply a suitable regularization (deformation) to each of these formal expressions 
so that the formal expression can be viewed as a singular value of the deformation function. 
Expanding around the singular point gives a Laurent series in k[£ -1 ,£|]. On the other hand, 
the algebraic combinatorial structure of the formal expressions, inherited by the deformation 
functions, can be abstracted to a free object in a suitable category. This free object parameterizes 
the deformation functions and often gives a Hopf algebra H. Thus the parametrization gives a 
morphism (p : H — » k[e _1 ,£]] in the suitable category. Upon applying the Algebraic Birkhoff 
Decomposition, we obtain (f>+ : H —> k[[e]] which, composed with e i-» 0, gives us well-defined 
values in k. 

4.2. Atkinson factorization. The following is the classical result of Atkinson. 

Theorem 4.2. (Atkinson Factorization) Let (R, P) be a Rota-Baxter algebra of weight A ^ 0. Let 
a G R. Assume that be and b r are solutions of the fixed point equations 

(34) b { = l+ P(b c a), b r = 1 + (id R - P)(ab r ). 
Then 

b c (l +Aa)b r = 1. 

Thus 

(35) 1 + Aa = b- c l b; 1 

ifbf and b r are invertible. 

We note that the factorization (l35l) depends on the existence of invertible solutions of Eq. (|34|) 
that we will address next. 

Definition 4.3. A filtered k-algebra is a k-algebra R together with a decreasing filtration R n , n > 
0, ofnonunitary subalgebras such that 

[J Rn = R, RnRm ^ R,i+m- 
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It immediately follows that Ro = R and each R n is an ideal of R. A filtered algebra is called 
complete if R is a complete metric space with respect to the metric defined by the subsets {R n }. 
Equivalently, a filtered k-algebra R with {R n } is complete if C\ n R n = and if the resulting embed- 
ding 

R^R:= limR/R n 

of R into the inverse limit is an isomorphism. 

A Rota-Baxter algebra (R, P) is called complete if there are submodules R n c R,n > 0, such 
that (R,R n ) is a complete algebra and P(R n ) Q R n . 

Theorem 4.4. (Existence and uniqueness of the Atkinson factorization) Let (R, P, R n ) be a com- 
plete Rota-Baxter algebra. Let a be in R\. 

(a) The equations in ffifli) have unique solutions b c and b r . Further b c and b r are invertible. 
Hence Atkinson Factorization 071) exists. 

(b) If A has no non-zero divisors in R\ and P 2 = —AP (in particular if P 2 = —AP on R), then 
there are unique C{ e 1 + P(R) and c r e 1 + (id# - P)(R) such that 

1 + Aa = c ( c r . 



4.3. From Atkinson factorization to algebraic Birkhoff decomposition. We now derive the 
Algebraic Birkhoff Decomposition of Connes and Kreimer in Theorem 14.11 from Atkinson Fac- 
torization in Theorem 14.41 Adapting the notations in Theorem 14.11 let H be a connected filtered 
Hopf algebra and let (A, Q) be a commutative Rota-Baxter algebra of weight A = -1 with Q 2 = Q, 
such as the pair (A, Q) in Theorem 14. II (see also Example 12 .41 ) . The increasing filtration on H in- 
duces a decreasing filtration R n = {/ e Hom(//,A) | f(H"~ l ) = 0},n > on R := Kom(H,A), 
making it a complete algebra. Further define 

P-.R^R, P(f)(x) = Q(f(x)), f e Hom(H, A), x 6 H. 

Then it is easily checked that P is a Rota-Baxter operator of weight -1 and P 2 = P. Thus 
(R,R n , P) is a complete Rota-Baxter algebra. 

Now let (p : H — » A be a character (that is, an algebra homomorphism). Consider e—<f> : H — » A. 
Then 

(e - = eOtf) - 0(1//) = 1// - 1// = 0. 

Thus e - (p is in R\ . Take e - to be our a in Theorem l4.41 we see that there are unique q g /'(^i) 
and c r e P(Ri) such that 

= c e c r . 

Further, by Theorem 14.21 for b c = c~ x , b e = e + P{b e ★ (e - <p)). Thus for x £ kere = kere, we 
have 

ZMX> = P(b { ★ (e - cp))(x) 

= Q^ a W^ e ~ ^)( a (2))) 

W 

= G(Wljr)(c - 0)W + J] ^(Jt'Xc - 0)(Jc") + - 0)(l ff ) 

(a) 

= -G(0(*) + £&K*W'))- 
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In the last equation we have used e(a) = 0, e{a") = by definition. Since bc(\n) = 1//, we see 
that b e = 0_ in Eq. (1321) . 
Further, we have 

c r = c jV = M = -(/>) + b c = -b c (e -</>) + e + P(b c (e -</>)) = e - (id - P)(b ( (e - <p)). 
With the same computation as for b ( above, we see that c r = 0+ in Eq. (1331) . 

5. Heat-kernel type regularization approach to the renormalization of MZVs 

To extend the double shuffle relations to MZVs with non-positive arguments, we have to make 
sense of the divergent sums defining these MZVs. For this purpose, we adapt the renormalization 
method from quantum field theory in the algebraic framework of Connes-Kreimer recalled in 
the last section. We will give three approaches including the approach in this section using a 
heat-kernel type regularization, named after a similar process in physics. Since examples and 
motivations of this approach can already be found elsewhere lT3T)ll3~6l I3TTI , we will be quite sketchy 
in this section. More details will be given to the two other approaches in Section [6] 

5.1. Renormalization of MZVs. Consider the abelian semigroup 

(36) 9Jf = {p|(s,r)eZxR >0 } 

with the multiplication 

[ r ] ' [/] — [ r + / ]■ 
With the notation in Section [2721 we define the Hopf algebra 

h w := MSc,i(cart) 

with the quasi-shuffle product * and the deconcatenation coproduct A in Section [2721 For w, = 
[ '] 6 9JI, i = 1, • • • , k, we use the notations 

'i 

w = (w u ...,w k ) = [*'"'• '*"] = U, where s = (s u . . . , s k ), ?= (r u ...,r k ). 

' 1 j ■ ■ ■ •> ' k ' 

For w = [t] e W k and e e C with Re(e) < 0, define the directional regularized MZV: 



(37) Z([%s)= J] 



S\ St 

n, ■ ■ ■ n. 

It converges for any [J and is regarded as the regularization of the formal MZV 

(38) m= y 1 „ 

ni>-»jjt>0 1 fc 

which converges only when Sj > and si > 1. It is related to the multiple polylogarithm 

z" 1 ■ ■ ■ z k 

U Sl _ Sk ( Zl ,...,z k ) = J] n Si ---n Sk 

by a change of variables Zi = e ns , 1 < i < k. 

This regularization defines an algebra homomorphism 11361 : 

(39) Z:H m ^C[T][[e,e- l l 
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In the same way, for 

(40) ST ={[ s r ] |(5,r)GZ< xR >() }, 

Z restricts to an algebra homomorphism 

(41) Z:fK 9 j-->2?:=C[[e,e _1 ]. 

Since both (C[r][e _1 , e]], IT) and (C[e _1 ,e]],IT), with II defined in Example 12.41 are commu- 
tative Rota-Baxter algebras with IT 2 = IT, we have the decomposition 

z = z: 1 ★ z + 

by the algebraic Birkhoff decomposition in Theorem 14. H and obtain 

Theorem 5.1. ( Il36l [371 ) The map Z + : "K m — > C|T][[£]] is an algebra homomorphism which 
restricts to an algebra homomorphism Z + : %m- — > C[[e]]. 

Because of Theorem (5TTJ the following definition is valid. 

Definition 5.2. For s = (s u . . . , s k ) e Z k and r = (r l5 . . . ,r k ) e R k >0 , define the renormalized 
directional MZV by 

(42) <r($) = UmZ + ($;s). 

' £— »0 

Here ris called the direction vector. 

As a consequence of Theorem (5TTJ we have 
Corollary 5.3. The renormalized directional MZVs satisfy the quasi-shuffle relation 

(43) CdfMiU) = * £])■ 
Here the right hand side is defined in the same way as in Eq. ((#]). 
Definition 5.4. For s e Z k >0 U Z* , define 

( 44 ) M = ^(i[/ +S ^ 

where, for s = (si, ■ ■ ■ ,Sk) and 5 e R >0 , we denote \s\ = (\s\\, • • • , \sk\) and \s] + 6 = (\s\\ + 
5, • • ■ , \Sk\ + 6). These £(s) are called the renormalized MZVs of the multiple zeta function 
£(ui, ■ ■ ■ ,uu) at s. 

Theorem 5.5. (35| 

(a) The limit in Eq. ft44\) exists for any s*= (s\, ■ ■ ■ , Sk) 6 Z k >0 U Z k <0 . 

(b) When s { are all positive with S\ > 1, we have £([^]) = £(s) independent ofre Z k >Q . In 
particular, we have %(s*) = £(s). 

(c) When s t are all positive, we have £(s) = Further, l(s) agrees with the regularized 
MZV Z*JT) defined by Ihara-Kaneko-Zagier [|46l . 

(d) When Sj are all negative, we have g(s*) = £([_ J) = lim £(["!>]). Further, these values are 
rational numbers. 



DOUBLE SHUFFLE RELATIONS AND RENORMALIZATION OF MULTIPLE ZETA VALUES 21 

(e) The value <f(j) agrees with £(s*) whenever the latter is defined by analytic continuation. 
Furthermore, 

(f) the set (l{s)\s e Z* } satisfies the quasi-shuffle relation; 

(g) the set (l{s)\s e Z k <0 } satisfies the quasi-shuffle relation. 

The following table lists £(-a, -b) for < a, b < 6. 
(45) 



<f (-a, -fc) 


a = 


a = 1 


a = 2 


a = 3 


a = 4 


a = 5 


a = 6 


- n 


3 
8 


l 

12 


l 

120 


l 

120 


l 

252 


l 

252 


l 

240 


h - 1 


1 

24 


1 

288 


1 

240 


83 


1 

504 


3925 


1 

480 


64512 


2239488 


b - 2 


1 

120 


1 

240 


o 


1 

504 


319 


1 

480 


2494519 


437400 


1362493440 


b = 3 


1 

240 


71 


1 

504 


1 


1 

480 


114139507 


1 

264 


35840 


28800 


139519328256 


b = 4 


1 

252 


1 

504 


319 


1 

480 





1 

264 


41796929201 


437400 


26873437500000 


b = 5 


1 

504 


32659 


1 

480 


21991341 


1 

264 


1 


691 


15676416 


25836912640 


127008 


65520 


b = 6 


1 

240 


1 

480 


2494519 


1 

264 


41796929201 


691 





1362493440 


26873437500000 


65520 



5.2. The differential structure. The shuffle relation for convergent MZVs from their integral 
representations does not directly generalize to renormalized MZVs due to the lack of a suitable 
integral representation. However a differential variation of the shuffle relation might exist for 
renormalized MZVs. One evidence is the following differential version of the algebraic Birkhoff 
decomposition 11371 for renormalized MZVs and further progress will be discussed in a paper 
under preparation. We first recall some concepts. 

(a) A differential algebra is a pair (A, d) where A is an algebra and d is a differential op- 
erator, that is, such that d{xy) = d(x)y + xdiy) for all x,y e A. A differential algebra 
homomorphism / : (A\, d\) — > (A 2 ,d 2 ) between two differential algebras (A\,d\) and 
(A2, di) is an algebra homomorphism / : A\ — » A 2 such that / o d\ = d 2 /. 

(b) A differential Hopf algebra is a pair (H, d) where H is a Hopf algebra and d : H —> H is 
a differential operator such that 

(46) A(d(x)) = J] (d(x {i) ) (g) x (2) + x (1) (g) d(x {2) )). 
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(c) A differential Rota-Baxter algebra is a triple (A, IT, d) where (A, II) is a Rota-Baxter 
algebra and d : R — » i? is a differential operator such that P o d = d o P. 

Theorem 5.6. (Differential Algebraic Birkhoff Decomposition) 0711 Under the same assump- 
tion as in Theorem \4.1\ if in addition (H, d) is a differential Hopf algebra, (A, II, d) is a commu- 
tative differential Rota-Baxter algebra, and (p : H — > A is a differential algebra homomorphism, 
then the maps (j>- and (p + in Theorem \4.1\ are also differential algebra homomorphisms. 



Theorem 5.7. ( 071 ) 

(a) For [ s ] £ 9K, define d([ s r ] = r[ s ~ ]. Extend d to "K m = ® k>Q (kWf k by defining, for 
a:=ai®---®a k £(km) IBk , 



Then (IKU, d) is a differential Hopf algebra. 

(b) The triple (C[e _1 , e]], II, 4-) is a commutative differential Rota-Baxter algebra. 

(c) The map Z : -> C[[e, e" 1 ] defined in Eq. d?7J) is a differential algebra homomorphism. 

(d) The algebra homomorphism Z + : "K^- - > C[[e]] in Theorem \5.1\ is a differential algebra 
homomorphism. 

6. Renormalization of multiple zeta values seen as nested sums of symbols 

We present two more approaches to renormalize multiple zeta functions at non-positive inte- 
gers, both of which lead to MZVs which obey stuffle relations. Like the renormalization method 
described in the previous section, they both give rise to rational multiple zeta values at non- 
positive integers and we check that the two methods yield the same double multiple zeta values at 
non-positive integer arguments. This presentation is based on joint work of one of the authors with 
D. Manchon 051 in which multiple zeta functions are viewed as particular instances of nested 
sums of symbols and where the algebraic Birkhoff decomposition approach is used to renormal- 
ize multiple zeta functions at poles. Here, we furthermore present an alternative renormalization 
method based on generalized evaluators used in physics ll67ll . 

6.1. A class of symbols. For a complex number b, a smooth function / : R - {0} — » C is called 
positively homogeneous of degree b if = for all t > and £ e R. 

The symbols which were originally defined on R" are now defined on R which is sufficient for 
our needs in this paper. We call a smooth function <x : R — > C a symbol if there is a real number 
a such that for any non-negative integer y, there is a positive constant C y with 



For a complex number a and a non-negative integer j, let cr a -j : R - {0} — » C be a smooth and 
positively homogeneous function of degree a- j. We write <x ~ Y*J=o °~a-j if, for any non-negative 
integer N and non-negative integer y, there is a positive constant C y>N such that 



(47) 




|<9V(£)| < C y {\ + m>-\ V£ e R. 




& o-(Z)-Y,o- a -tf) <C 7tN (l + \$) 



Re(a)-iV-l-7 



V^R-(O), 



V ;'=0 ) 

where Re(<2) stands for the real part of a. 
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For any complex number a and any non-negative integer k, a symbol cr : R — » C is called a 
log-polyhomogeneous of log-type k and order a if 

k oo 

(48) oig) = J] log' |£|, o-tf) ~ Yj 

with cr a _ ; / (^) positively homogeneous of degree a - j. 

Let §>"' k denote the linear space over C of log-polyhomogeneous symbols on R of log-type k 
and order a. Then we have §>"- k c § a < k+l . Let §>*' k denote the linear span over C of all §> a,k for a e C. 
Then S*'° corresponds to the algebra of classical symbols on R. We also define 



*,k 

o 

k=0 



which is an algebra for the ordinary product of functions filtered by the log-type 11521 since 
the product of two symbols of log-types k and k! respectively is of log-type k + k'. The union 
U«€Z UT=o & a ' k 1S a subalgebra of §*•*, and [j aeZ S a, ° is a subalgebra of S*'°. 

Let f S >a ' k be the algebra of positively supported symbols, i.e. symbols in § a,k with support in 
(0, +oo) so that they are non-zero only at positive arguments. We keep mutatis mutandis the above 
notations; in particular is a subalgebra of the filtered algebra T* *. 

For cr 6 7 a ' k we call fp <r(£) := cr Q; o(£) the finite part at zero (so named since it it reminiscent 

of Hadamard's finite parts) of such a symbol cr which corresponds to the constant term in the 
expansion. 

The following rather elementary statement is our main motivation here for introducing log- 
polyhomogeneous symbols. 

Proposition 6.1. Il54ll The operator I defined in ([2]) on the algebra C[0, oo) by 

f^^i(fW = ^f(t)dt ) j 

maps y*-^ 1 to "P** for any positive integer k. 

By Proposition 16.11 for any cr in 7*' k , the primitive 7(cr)(^) has an asymptotic behavior as 
£ — > oo of the type (|48l ) with k replaced by k + 1 . The constant term defines the cut-off regularized 
integral (see e.g. Il52l ): 

-I cr{t)dt := fp I cr(t)dt. 

JO £->oo Jo 

6.2. Nested sums of symbols and their pole structures. 

6.2.1. Nested sums. Recall that the operator / on 3 *'* defined by Eq. (0 satisfies the weight 
zero Rota-Baxter relation (OQ). On the other hand the operator P defined by Eq. © satisfies the 
Rota-Baxter relation with weight A = -1 and the operator Q = P - Id in Eq. © satisfies the 
Rota-Baxter relation with weight A = 1 . 

The Rota-Baxter operators P and / relate by means of the Euler-MacLaurin formula which 
compares discrete sums with integrals. For cr e 7*'* the Euler-MacLaurin formula (see e.g. IT3810 
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reads: 



P(cr)(N) - I(cr)(N) = -cr(N) + }]^ k ' l) (N) 



k=2 



(49) + 1 f B 2K+l (x)o- (ZK+1 \x)dx. 

(IK + I)! Jo 

with B k {x) = B k {x- [x]). Here B k (x) = 2?=o (;) Bk-i x k are the Bernoulli polynomials of degree 
k, the B t being the Bernoulli numbers, defined by the generating series: 

_L_ = y 5^. 

e'-l ^ i\ 

Since B k (\) = B k for any k > 2, setting jc = 1 we have 

(50) 5, = ^ ( k \ B k _i = £ (*] P„ V£ > 2. 

i=o \ 1 ' 1=0 ^ Z '' 

The Euler-MacLaurin formula therefore provides an interpolation of P(cr) by a symbol. 

Proposition 6.2. Il55l For anv a e D 5 "'*, ?/ze discrete sum P(cr) can be interpolated by a symbol 
P(o-) in y« +u+1 + 9 0M1 (i.e. P(o-)(n) = P(a)(n) = 2X0 a(k), "in e N) such that 

P(o-) - 1(a) e T' k . 

The operator Q :=P-Id: T' k T ,+lM1 + T Ml interpolates Q. 

By Proposition 16.21 given a symbol o~ in 7"' k , the interpolating symbol P(o~) lies in y a+u+1 + 
rpo,£+i j t follows that the discrete sum P(cr)(N) = P(cr)(N) has an asymptotic behavior for large 
N given by finite linear combinations of expressions of the type (|48T) with k replaced by k + 1 and 
a by a + 1 or 0. Picking the finite part, for any a e 7*'* we define the following cut-off sum: 

oo N 

(51) ¥>:= fp P(a)(N) = fp Ya(k), 



which extends the ordinary discrete sum Xo° on L l -symbols. If <x has non-integer order, we have 

OO CO 

£ °" = fp Zf=cf cK^) f° r an Y integer ^, so that in particular £ <x = fp Q(cr)(N) since the 

iV-»oo iV^oo 

operators P and 2 only differ by an integer in the upper bound of the sum. 

With the help of the interpolation map described in Proposition I6.2[ we can assign to a tensor 
product a := a\ <S> ■ ■ ■ ®a k of (positively supported) classical symbols, two log-polyhomogeneous 
symbols defined inductively in the degree k of the tensor product, which interpolate the nested 
iterated sum 



2^ o-\(n{) ■ • ■ cr k (n k ) = cr, P( • • • cr fc _ 2 P(o- fe -i P(a k ))..)j, 
2^ o-i(n x ) ■ ■ ■ cr k (n k ) = ctj g( • • • cr^ <2(cnt-i P(cnO)...). 

0<«t <n,t- 1 < • ■ ■ <«2 <« 1 

In the following we will only consider the second class of symbols, including their regularization, 
renormalization and application to multiple zeta values. A parallel approach applies to the first 
class of symbols with application to non-strict multiple zeta values in Eq. (1221) |j6TH75l . 
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Theorem 6.3. 11551 Given o~j e J"*" , i = 1, . . .,k, setting cr := cri <g> • • • <8) cry the function cr 
defined by: 

(52) 5= := o-i fi( • • ■ o-*_ 2 G(tr t ))...) 

which interpolates nested sums in the following way: 

cr{n x ) = 2j cr 1 (n 1 )---cr fc (n fc ), Vn x g N, 

0<«t<fi,t_i<---<n2<ni 

/zes in T*' 4-1 a* linear combinations of (positively supported) symbols in , y> a i+ -+ a j+j- l J- l t j £ 
{l,...,k}. 

On the grounds of this result, we define the cut-off nested discrete sum of a tensor product of 
(positively supported) classical symbols. 

Definition 6.4. For o~i, . . . , cr k e T* ,Q and cr := o~\ ® ■ ■ ■ <8» cr k we call 

Chen oo 

J]o-: = J]o : = cr^m) • • • o- k (n k ) 

< 0<%<— oil 

the cut-off nested sum of / = cr\ ® • • • ® cr k . 

6.2.2. 77ze po/e structure of nested sums of symbols. To build meromorphic extensions, we com- 

oo 

bine the cut-off sum £ introduced in (|5T| ) with holomorphic deformations of the symbol in the 

o 

integrand. 

A family {a(z)} z en in a topological vector space A which is parameterized by a complex domain 
CI, is holomorphic at z g Q if the corresponding function f : CI —> A admits a Taylor expansion 
in a neighborhood N zo of z 

(z-zof 



a(z) = J]a ik \zo) 



kl 

which is convergent, uniformly on compact subsets of N zo (i.e. locally uniformly), with respect 
to the topology on A. The vector spaces of functions we consider here are C(R, C) and C°°(R, C) 
equipped with their usual topologies, namely uniform convergence on compact subsets, and uni- 
form convergence of all derivatives on compact subsets respectively. 

Definition 6.5. Let k be a non-negative integer, and let Q be a domain in C. A simple holomor- 
phic family of log-polyhomogeneous symbols cr(z) g §*' k parameterized by CI is a holomorphic 
family <x(z)(£) := a(z, £) of smooth functions on R such that: 

(a) the order a : CI — » C is holomorphic on CI, 

(b) oizW = llo<ri(& log' £| with 

<r,(z)(0 ~ J] ^Uw,/(£>- 

Here cr^)^)^ positively homogeneous of degree a(z) - j, 

(c) for any positive integer N there is some positive integer K N such that the remainder term 

k K N 

<r w {zm ■= oizXO "EE a ^HM) log' If I = I"") 

/=0 ;=0 
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is holomorphic in z e Q as a function of £ and verifies for any e > the following 
estimates: 

<ffi<r {N) (z)(i;) = o(\Z\- N -^) 

locally uniformly in z e Q for k e N and fi eW. 

A holomorphic family of log-polyhomogeneous symbols is a finite linear combination (over C) 
of simple holomorphic families. 

It follows from the Euler-MacLaurin formula (see e.g. Il38l [551 ) that for any holomorphic 
family cr(z) of symbols in "?*•*, we have 

V cr(z)(n) = A oizW d£ + C(cr(z)) 

with z i-> C(cr(z)) a holomorphic function at zero. Hence, z 2 cr (^)( n ) an d z i-> -t cr(z)(£) d£ 



have the same pole structure. Results by Kontsevich and Vishik ll49l for classical symbols and 
their generalization by Lesch ll52l to log-polyhomogeneous symbols, and relative to the pole 
structure of cut-off integrals of holomorphic families of symbols, therefore carry out to discrete 
cut-off sums of holomorphic (positively supported) log-polyhomogeneous symbols. Let us briefly 
recall the notion of holomorphic regularization inspired by 113911 . 

Definition 6.6. A holomorphic regularization procedure on §* * is a map 

ft : S*'* -> Hol n (§*'*) 

/ >-» {o-{z) = o- f (z)} ze n 

where Q. is an open subset of C containing 0, and Holn (§*'*) is the algebra of holomorphic families 
in §* * , such that for any / 6 §*•*, 

(a) o-(0) = /, 

(b) the holomorphic family <r(z) can be written as a linear combination of simple ones: 

k 

o-(z) = ^ o-j(z), 

M 

the holomorphic order aj(z) of which verifies Re(a'j(z)) < for any z e Q. and any 
J'e{l, ...,£}• 

A holomorphic regularization ft is simple if, for any log-polyhomogeneous symbol cr e $ a ' k , the 
holomorphic family ft(cr) is simple. Since we only consider simple holomorphic regularizations, 
we drop the explicit mention of simplicity. 

A similar definition holds with suitable subalgebras of §*•*, e.g. classical symbols §*'° instead 
of log-polyhomogeneous. Holomorphic regularization procedures naturally arise in physics: 

Example 6.7. Let z i-> t{z) e §*'° be a holomorphic family of classical symbols such that r(0) = 1 
and t(z) has holomorphic order a(z) with Re(a'(z)) < 0. Then 

ft : cr h-> cr(z) := ctt(z) 

yields a holomorphic regularization on §* * as well as on §*'°. Choosing r(z)(^) := x(%) + (1 _ 
X(£))(H(z) |^| _z ) where H is a scalar valued holomorphic map such that H(0) = 1, and where x 
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is a smooth cut-off function which is identically one outside the unit interval and zero in a small 
neighborhood of zero, we get 

B(<r)(zX0 = *(£M£> + (1 -x(&)(H(z) oi& {ZD. 
Dimensional regularization commonly used in physics is of this type, where H is expressed in 
terms of Gamma functions which account for a "complexified" volume of the unit sphere. When 
H = 1, such a regularization CR is called Riesz regularization. 

Proposition 6.8. Given a holomorphic regularization : o~ i-» o~(z) on CP*'*, for any o~ e 

oo 

CP*'* the map z i-» £ o~(z) w meromorphic with poles of order < k + \ in the discrete set 
o 

a" 1 ({-1,0, 1, 2, • • • }) whenever cr(z) z's a holomorphic family with order a(z) such thatRe.(a'(z)) + 
Ofor any z in Q. 

Let Q. c C be an open neighborhood of 0. Given symbols cr\, • • • ,crt e CP*' , and a holo- 
morphic regularization CR which sends <r, to cr,(z) with order ar ; (z), z £ Q, we build holomorphic 
perturbations in the complex multivariable z := (zi,- ■ • ,Zk) € & k of the symbols cr introduced in 
<E2]>: 

<r(z) := o-i(zi)fi(-- -o-k-zizk-z) Q(o- k -i(z k -i) Q(o- k (z k )))...y 

By Theorem |6.3[ these are linear combinations of log-polyhomogeneous symbols of log-type 
j — 1 and order ari(zi) + • • • + ar/(z/) + y - 1, y e {1, . . . , k}. Applying Proposition [6]8] to each of 
these symbols provides information on the pole structure of nested sums of (positively supported) 
classical symbols reminiscent of the pole structure of multiple zeta functions lfl"l [2~5H73ll . 

Theorem 6.9. Fix symbols o~\, ■ ■ ■ ,crk e CP*' and a holomorphic regularization CR which sends 
o~i to o~i(z) with order ar,-(z). 

(a) The map 

Chen 

(Z\,---,Zk) ^CTi(Zi)® •• •®£Tjfc(Zjfc) 

< 

Z5 meromorphic with poles on a countable number of hypersurfaces 

i 

Yj anfa) e -j + N , 
i=\ 

of multiplicity j varying in {1, • • • , /c}. //ere No stands for the set of non-negative integers. 

(b) Le? cr(z) := cr^z) <8> • • • <8> o"*(z) wzY/i z 6 Q. Assume that the orders ar;(z) of the o~i's are 
nonconstant affine with a'-(Qt) = —qfor any j in {1, • • • , k] and some positive real number 

Chen 

q. The map ^ <x(z) is meromorphic on Q. with poles z e (Z/=i <*i(0) + J ' ~ No)/(g /) 

< 

of order < j. 

(c) 7/'Re(Q'i(zi) + • • • + (Xj(Zj)) < -j for any j e {1, . . . ,k}, the nested sums converge and boil 
down to ordinary nested sums ( independently of the perturbation). Setting cr = cr\®- ■ -®crk 
we have: 

Chen, 3? Chen Chen 

< < < 
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6.3. A twisted holomorphic regularization. We now take A to be a subalgebra of T* ° equipped 
with the ordinary product on functions. Any holomorphic regularization 31 on A with parameter 
space QcC induces one on the tensor algebra T(A): 

^(o-i ® • • • ® o- k )( Zl , ■■■ ,z k ):= ft(o-i)(zi) ® • • • ® ^(o- ft )fe). 

It is compatible with the shuffle product 

31 ((o"! ® • • • ® cr*) in (o- M ® • • • ® cr^,)) = IR (a"! (8) • • • ® a k ) mJi (cr A . +1 ® • • • <g> cr k+ i) 

for any o- t e A, i e {I, ■ ■ ■ ,k + I}. 

Remark 6.10. Note that 3i(cri incr 2 Xzi,z 2 ) + ^(o"i)(zi) m3i(a 2 )(z 2 ) whereas 3i(cri mcr 2 )(zi,z 2 ) = 
(R(o- x )wR{o- 2 ))(zuZ 2 ). 

Let 

8 k : C -> C^, z i-> z • l 8 *, 
be the diagonal map 5 : C h T(C) and 5* the induced map on tensor products of holomorphic 
symbols 

51 : T (Holn (A)) -» Hol n (7(A)) , a h-> cr o 
The regularization induces a one parameter holomorphic regularization: 

(6* o 3t) (o-! ® • ■ • ® trjOCz) = )(z) ® • • • ® #(cr A )(z) 
compatible with the shuffle product: 

(<T o ft) ((o-j ® ••• ® 0-^) m (cr fc+1 ® • • • ® o- k+l )) 

(53) = (<T o ft) (en ® • • • ® a*) in (5* o ft) (cr, +1 ® • • • ® cr k+l ) , 

for any cr, ■ e A, i 6 { 1 , • • • , k + I}. 

Twisting it by Hoffman's isomorphism in Theorem 12.81 yields a holomorphic regularization 
(s* o ft) * (denoted by ft* in (551) on r(A): 

(5* o ft) * : = exp o (tf* o ft) o log, 
which is compatible with the stuffle product: 

(54) (<T o ft)* (cr * r) = (cf o ft)* (cr) * (<f o ft)* (r), Vcr, t £ T (A) . 
Consequently, the following regularization 

(55) T(o-i ® • • • ® cr k )(zu ■ ■ ■ ,z k ) = exp oft o log^ ® • • • ® cr k ) (zu ■ ■ • ,z*) 
is compatible with stuffle relations after symmetrization in the complex variables zt 

(56) (ft*(cr*r)) = (ft*(cr) * ft*(r)) , Vo-,t£T(A), 

V /sym V /sym 

where the subscript sym stands for symmetrization 

/sym(Zl> - " - ,Z k ) '■= — ^ /(Z T (1), • • • , Z T (fc)), 

over all the complex variables z\, ■ • • , z k +i if cr is a tensor of degree k and r a tensor of degree /. 
Setting Zi = • • • = Zjt+? = Z in (l56l) yields back (l54l) so that (l56l) can be seen as a polarization of 
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(|54|) . Given symbols a\, ■ • ■ , cr k in CP*' , and a holomorphic regularization CR, : cr 1— » <x(z), sending 
<Ji to (Ti(z) with order we are now ready to build a map 

Chen 

(zi, ■•• 2 ® • • • ® o*)(zi, • • • .**)» 

< 

which, by Theorem 16.91 is meromorphic with poles on a countable number of hypersurfaces 

J] e -j + No, 

1=1 

with multiplicity j varying in {1, • • • , k}. In particular, if the holomorphic regularization sends 
a symbol <x to a symbol cr(z) with order a(z) = a(0) — qz for some positive real number q, the 
hypersurfaces of poles are given by 

so that hyperplanes of poles containing the origin correspond to YJ = \ Z; = each of which with 
multiplicity j varying in {1, • • • , k}. 

6.4. Meromorphic nested sums of symbols. Let Mer (C) denote the germ of meromorphic 
functions in a neighborhood of zero in the complex plane and let Holo(C) be the germ of holo- 
morphic functions at zero. We consider the (Grothendieck closure of the) tensor algebra 

r(Mer (C)) = Co^'(Mer (C)) 

over Mer (C) and its subalgebra r(Hol (C)) := ©£l o r*(Holo(C)) where we have set r*(Mer (C)) := 
®*Mer (C), r*(Hol (C)) := ®*Hol (C), and where <g> stands for the Grothendieck closure. They 
come equipped with the product: 

C/i ® • • ■ ®/*) ® • • •»/*+/) = /i ® • • • ® ■ ■ • ®/*+/- 

We consider the following linear extension of 7*(Mero(C)) which corresponds to germs at zero of 
meromorphic maps in severable variables with linear poles. Let £Mero(C°°) := ©j^ZLMeroCC*) 
where 

{m 
lJfi° L i\ ft e Mer (C), L,-e(c*) 

or equivalently, 

£Mer (C*) := L, ■ ■ ■ ,z k ) h» K % ' ' ' Zk) - 

Setting m = k and L,(zi, • • • , Zk) = Z\ yields a canonical injection 



/i E Hol (<C*) , m L e 



i : r(Mer (C)) -> £Mer (C) 



(zi,--- ,z*) >-» l^/io^Kzi, -- ,z*) 
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and the tensor product on r(Mer (C)) extends to £Mer (C°°), by 



(57) 



i=\ 



(Z\, ■■■ ,Zk) ^> Y\fi ° L i(Z\,-- - >Zk) • (ZU- ■ ■ ,Zl) ^ J~[ fm+j ° L i+j(Zl, ■ ' 1 ,Zl) 



(Z\, ■■■ ,Zk,'" >Zk+l) >-> Pj/i ° Lj(Zi,- ■■ ,Zk) \ \fm+j ° Lm+j(Zk+W ' ' , 
'=1 7=1 

which makes it a graded algebra. 

Specializing to linear forms L k := {L e (C fe J | 37 c {1, • • • ,k], L(zi, ■ ■ ■ ,Zt) = JjjejZj] 
gives rise to a subalgebra £M (C°°) := ®^ =l LM Q (C k ) c £Mer (C°°) defined by: 

h(zi, ■ ■ ■ ,z k ) 



£M () (C) :=Uzu--- ,z k )» 



For future use, we consider the map 6* : £Mo(C*) — > Mer (C) defined by 

(5* : £Mo(C fc ) -» Mer (C), f»fo5 k , 



h e Hol (<C*) , m L 6 



induced by the diagonal map 5 : C i-» r(C) previously defined. 



Chen 



By definition of the twisted regularization CR*, the expressions £ 



Chen 



are linear combinations of expressions of the type ■£ Ti(«i) <S> • • • <8> Ti(uk) with symbols Tj(uj) 

< 

built from products of the cr^zO's. It therefore follows from Theorem 16.91 that the functions 

Chen _ 

(zi, • • • ,Zk) £ 5i*(cr\ ® ■ ■ ■ <g> cr^Xz!, • • • ,z fe ) lie in £Mo(C*). Since the stuffle relations are 

< 

satisfied for convergent nested sums, given two tensor products <x = <T\ ® • • • <g> cr k and r = 
• • <g> Ti of symbols in CP, setting cr,-(z,) := CR(<x ; )(z ( ), for Re(z,) sufficiently large we have: 

Chen 

£ (W (o-i ® • • • ® CTfc) ( Zl , • • • , Zjfc )) * (ft* (n ® ■ • ■ ® T,) (Z t+ i, • • • ,z k+ d) 



< 

/Chen 



£ft*( 



crjtXzi, • • • ,Z*) 



/Chen 



T/XZfe+i* - - - ,Zt+/) 



V < / V < 

By analytic continuation (see for example [|29l . in particular the Identity Theorem in Chapter 1, 
Section A, or 05)), this holds as an identity of meromorphic functions. Since 



(ft* ((o-i 9 ■ ■ ■ <8 o- h ) * (n ® ■ ■ ■ ® t,)) (zi, • • • ,Z M )) 



sym 



= ((ft*(0"l <g> • • • ® CTfe) * ft* (Tl ® • • • ® T)) (Zfa-i , • ' • , Z k+ l)) s ■ 



symmetrization in the variables n yields 



/Chen 



(&Vi ® • ' • ® O-j) * ft*(Ti ® • • • ® T,)) (Zl, • • • ,Z k+ l) 



\ < 
/Chen 



sym 



/Chen 



£ ft*(cri (8) • • • ® cr^Xzj, • • • ,Z k ) £ ft*(Ti ® • • • ®T/XZk+l, • • • ,Z k +l) 



J \ < 



sym 
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This can be reformulated as follows. 

Theorem 6.11. Il55l Let A be a subalgebra of !P*'° and let R be a holomorphic regularization 
which sends a symbol cr to a symbol o~(z) with order a(z) = a(0) - qzfor some positive real 
number q. 

(a) The map 

T^ : (T (A) , *) -> (£Mo(C°°), •) 



/ Chen 

CT\ <8> • • • <8> CT k h-» 



(zi,--- ,z k )^ J] 91* (ai ®---®o- k )(zi,-- - ,Zk)\, 



satisfies the following relation: 

(^(o- * t)) = (T K (o-) • ¥ k (t)) 

V /svm V /s 



' sym v / sym 

which holds as an equality of meromorphic functions in several variables. Here, as before 
the subscript sym stands for the symmetrization in the complex variables Zu 
(b) After composition with 8* this in turn gives rise to a map 



ijr*- : (T (A) , *) -> Mer 

( Chen \ 

zt-»iTo£ ^*(o-i ® • • • ® cr fc )(z) 



(58) o"! (8) • • • (8) a- A i-> 



which is an algebra morphism. In other words, ip satisfies the relation: 

<AV * r) = <A V) • <A K (r), 
which holds as an equality of meromorphic functions in one variable. 

6.5. Renormalized nested sums of symbols. We want to extract finite parts from the mero- 
morphic functions in Theorem 16.111 while preserving the stuffle relations using a renormalization 
procedure. Renormalized evaluators inspired from generalized evaluators used in physics provide 
a first renormalization procedure. 

6.5.1. Renormalized nested sums via renormalized evaluators. We call regularized evaluator at 
zero on the germ Mer (C) of meromorphic functions around zero, any linear form on Mer (C) 
which extends the evaluation at zero ev : h i-» /i(0) on holomorphic germs at zero. The map evQ g 
defined by 

ev^ g :=ev o(/-n), 

where II : Mer (C) — > Mer (C) as defined in Example 12.41 corresponds to the projection onto 
the pole part of the Laurent expansion at zero, is such a regularized evaluator at zero. When we 
need to specify the complex variable z we also write ev^f . Following Speer [[671 we introduce 
renormalized evaluators which correspond to his generalized evaluators. 

Definition 6.12. A renormalized evaluator A on a graded subalgebra H = ® < ? Q r B k o/XMeroCC 00 ) = 
®™ =0 JLMer (C k ) equipped with the product • introduced in d57l) . is a character on B which is com- 
patible with the filtration induced by the grading and extends the ordinary evaluation at zero on 
holomorphic maps. Equivalently, 

(a) Compatibility with the filtration: Let H> K := ®f =0 ¥>k and Ak '■= A| K . Then Ak+i\ K = A%. 
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(b) It coincides with the evaluation map at zero on holomorphic maps: 

A lr(Ho l0( c» = eV 0" 

(c) It fulfills a multiplicativity property: 

A(/.g) = A(/)A(g), V/.geS. 
We call the evaluator symmetric if moreover for any f in H> k and r in we have 

Mfr) = A(/), Vr e I,, 
where we have set f T (zu ■■■ , Zi) ■= f(z T (i), ■■ , z T (k))- 

Example 6.13. Any regularized evaluator at zero A on Mer (C) uniquely extends to a renormal- 
ized evaluator A on the tensor algebra (T (Mer (C)) , <S>) defined by 

k 

i=\ 

Example 6.14. Any regularized evaluator A on Mer (C) extends to renormalized evaluators A 
and A' on £Mer (C°°) defined on LMzr {C k ) by 

A:=A ZI o---oA Zk , A' := A Zk o-..oA zl 

and to a symmetrized evaluator defined on LMer Q (C k ) by 

ASym: =^Z^ o ---°^ 

' T€L k 

where A Zj stands for the evaluator A implemented in the sole variable Zu the others being kept 
fixed. Their restrictions to T (Mero(C)) all coincide with A. 

Example 6.15. Take A := evj, eg , and set with the above notations 

ev™ := A; ev™' := A', ev™' sym := A sym , 

then given a holomorphic function /z(zi, Z2) in a neighborhood ofO and setting f(zi,Z2) '■= ^+'^ > 
we have 

ev- (f) = d in (0, 0); ev™'(/) = «0, 0); ev™^ (/) = g i W) + W), 0) = ev - § 6 * (/) , 
though in general, 

ren,sym . reg o* 

ev Q 3 ± ev 5 o 5 . 

Proposition 6.16. Let Abe a subalgebra of CP*' and let 01 be a holomorphic regularization which 
sends a symbol f to a symbol o~(z) with order a(z) = a(0) - qzfor some positive real number q. 
Let tbe a symmetrized renormalized evaluator on £Mo. The map 

¥ XE : (T (A) , *) -> C 

o"i ® • • • <g> o- k i — > £ o W^icri ® • • • ® o- k ) 

defines a character. In other words, satisfies the stuffie relation: 

v F K ' £ ((r * t) = ^ %t {(f) ■ 

Remark 6.17. Here, we use the fact that for a symmetrized evaluator A we have A(/) = A(/ sym ) 
where as before the subscript "sym" stands for the symmetrization in the complex variables zu 
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This proposition gives rise to renormalized nested sums of symbols 

Chen,K,£ 

o- 1 ®---®o- k := x ¥ xe (o- 1 ®---®cr k ) 



which obey stuffle relations: 



Chen,3?,£ 



/Chen,K,£ \ {Chen,%£ \ 



£ (^* r )= £ a £ r 



V < 



< / 



6.5.2. Renormalized nested sums via algebraic Birkhoff decomposition. On the other hand, the 
tensor algebra T{A) can be equipped with the deconcatenation coproduct: 



A(o-i 



cr k ) := Ef " 1 

7=0 



°>1 



which then inherits a structure of connected graded commutative Hopf algebra H2l . Using the 
convolution product ★ associated with the product and coproduct on T(A) and since Mer (C) 
embeds into the Rota-Baxter algebra C[e _1 , e]] we can implement an algebraic Birkhoff decom- 
position as in (f3TT) to the map if/ 31 in Eq. (1581) : 



associated with the minimal substraction scheme to build characters 

Proposition 6.18. Il55l Let Abe a subalgebra of 1 ?*' and let "Jibe a holomorphic regularization 
which sends a symbol a to a symbol o~(z) with order a(z) = a(0) - qz for some positive real 
number q. The map 

^ Bkk : (T(A),*) -> C 

o"! <g> • • • <g> o- k i-» i//^(0)(<Ti <g> • • • <g> o- k ) 

defines a character 

The map yields an alternative set of renormalized nested sums of symbols 

Chen,D?,Birk 

a- 1 ® • • • ® cr k := ^^(crj (8) • • • <g> cr fc ) 



which obey stuffle relations: 

Chen,X,Birk 

£ (<r*r)= E 



/"Chen,K,Birk "\ 
\ < ) 



rChen,K,Birk \ 
T 

/ 



E 



6.6. Renormalized (Hurwitz) multiple zeta values at non-positive integers. 
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6.6. 1 . An algebra of symbols. Since we consider both zeta and Hurwitz zeta functions, let us first 
observe that for any non-negative number v and any cr in y*' k , the map £ \-> t*cr(^) := <x(£ + v) 
defines a symbol in y*' k . 

Let A be the subalgebra of 7*-° generated by the continuous functions with support inside the 
interval (0, 1) and the set 

{cr 6 T'° 1 3v e [0, +oo), 3s e C, cr(£) = (f + vp when f > 1}. 

Consider the ideal N of A of continuous functions with support inside the interval (0, 1). The 
quotient algebra A = A/N is then generated by elements o~ s jV e with cr s v (^) = + v) _v for 
\£\ > 1 . For any v e R + the subspace A v of /I generated by {cr v v | 5 e C} is a subalgebra of A. We 
equip .A v with the following holomorphic regularization on an open neighborhood Q of in C: 

01: A v -> Hol n (A.) 

cr. s ,v ' * (Z i-> (1 cr. s , v +^ cr J+Z;V ) 

where ^ is any smooth cut-off function which is identically one outside the unit ball and vanishes 
in a small neighborhood of 0. 

Let W be the C- vector space freely spanned symbols indexed by sequences (wi, . . . , u k ) °f rea l 
numbers. In other words, W is T(W) where W = ® u€ ^Rx u where we identify x u with u for 
simplicity and set x u • x v = x u+v , u, v e R. We then define the stuffle product on W as usual in 
Eq. © or Eq. CLD with A = 1 . The map 

cr : W -> T(A V ), u = (u u ■ ■ ■ , u k ) i-> cr„ ;v := o- (B1> ... iBt;v) := cr Mi;v ® • ■ • <g> cr Mt;v 

induces a stuffle product on 7{A V ): 

0~u;v * 0~ u';v := °~ u*u';v 

As before, we twist the regularization 5L induced by H on T(A V ) by the Hoffman isomorphism 
(fT4l) to build a twisted holomorphic regularization 3i* in several variables which satisfies 

(£*(cr„; V ) * ft*(°V;v)) = (ft*(cW;v)) 
V /sym V /sym 

and a twisted holomorphic regularization c>* o 3?* in one variable compatible with the stuffle prod- 
uct: 

(6* o ft*(cr«;v)) * (5* o ft* (oV;v)) = 5* o ft*(cW;v). 

6.6.2. Multiple zeta values renormalized via renormalized evaluators. Let Q. be an open neigh- 
borhood of in C and let Ji : cr \-> {cr(z)} z en be the holomorphic regularization procedure on A 
previously introduced. The multiple Hurwitz zeta functions defined by: 

£(su . . . , Sk', v u ...,v k ) := ^((Tslv, ® • • • ® o- SktVk ) 

are meromorphic in all variables with poles^on a countable family of hyperplanes s\ + ■ ■ ■ + Sj ■ e 
] - 00, j] n Z, j varying from 1 to k. When vi = • • • = = v, we set 

£(su • • • , Sk', v) := £(su . . . , s k ; vi,...,v k ) 



When k — 2 and vi = ■ • ■ = v/ = v a more refined analysis actually shows that for some any negative real number 
v, poles actually only arise for s\ — -1 or si + S2 € {-2, —1,0, 2,4, 6, ■ • •}. 
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in which case they satisfy the following relations: 

U e (u*u';vj) = U E (u;v)£ & (u';v)) . 

V /sym V /sym 

The renormalized multiple Hurwitz zeta values derived from a symmetrized renormalized eval- 
uator £ on £Mo(C°°): 

i; E (s u ■ ■ ■ , s k ; vi,...,v*) := ¥ XE (o- SuVl <8> • • • <g> o- Sk , Vk ) 

denoted by ^^(tfi, . . ., Sk, v) when vi = • • • = V/t = v, satisfy stuffle relations in that case: 

£ E (u*u';v) = £ E (u;v)£ E (u';v). 

Let us compute renormalized values in the case k = 2 using a renormalized evaluator. For any 
a e R and m 6 N - {0} we introduce the notation: 

[a]j := a{a - 1) • • • (a - j + 1). 

We extend this to 7 = and j = -1 by setting: [a] := 1, := j^y. Combining Definition (1551) 

DT(cri <8>cr 2 )(zi,z 2 ) = crife) <g> cr 2 (z 2 ) - • cr 2 )(zi) + ^o-!(zi) • cr 2 (z 2 ) 

applied to the regularization 

R(o~i)(z)(x) = (x + v)~ Si ~ z of order = -s t - zu 

with the Euler-MacLaurin formula (1491) . and following 11551 (see the proof of Theorem 9), we 
compute 

£(Ji,J2;v)(zi,Z 2 ) = v P 3? (cr i ,, v «)cr i ,, v Xz 1 ,Z 2 ) 



Chen 

= ^0"l(Zl)®cr 2 (z2) + -cri(zi)cr 2 (z 2 ) - - (cr, o- 2 )(zi) 



< 

2./: 



^ 5; rp^— (fCsi + ft + zi + z 2 + 7 - 1; v) - ffa + Zu v)) 

1 1 

+ 5 2 + zi + z 2 ; v) - - £(si + s 2 + zi ; v) 



(27 2 + 1) 



f-^ {(n + vr Sl - zl B^~i(y)(y + vy 



Hence, for non-positive integers si = -a,\, s 2 = -a 2 and 2J 2 = a\ + a 2 + 2 we have: 

Z|/2 2 — Z2] j-1 
Bj (^(-ai - a 2 + zi + Zi + j - 1; v) - £(-a! + z x ; v)) 

7=0 J ' 

(59) +^(-«i - «2 + zi + z 2 ; v) - ^ ^"(-ai - a 2 + zi; v) 

[a 2 - z 2 L 2+2 
(fl 2 + 2)! , 

The last line on the r.h.s. is a holomorphic expression at zero on which all renormalized evaluators 
at zero vanish. The second line on the r.h.s is a linear combination of ordinary zeta functions at 



00 / r-n 

J]\(n + vr~ Zl ^ ~B a ~^{y)(y + vT 2 dy 
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negative integers which are holomorphic at zero. Any evaluator A at zero vanishes on these terms; 
indeed, we have A - a 2 + Z\ + z 2 , v) - £,{-ay - a 2 + zr, v)) = - a 2 , v) - £(-ai - 

a 2 ; v) = 0. Only when evaluated on the expression on the first line of the r.h.s can various 
evaluators differ. 

We want to implement the symmetrized evaluator at zero 



ren.sym 1 / reg reg reg 

2v n := - (ev j , _ A o ev_ _ n + ev_ _„ o ev 



• 2 V Z2=0 Zl=0 

introduced in Example 16.151 We first compute 

ev ^0 ( eV a=0 v )(Zl . Z 2 ))) 



=0 



reg \ 

Z2=0) 



= ev 



reg 
zi=0 



reg 



ev 



'«2 + l 



r ai+d2+2 

Z ^ 

V 7=0 



[«2 - Z2J7- 



- a 2 + zi +z 2 + 7 - 1; v) - £(-a, + Zi; v) 



= ev 



reg 
zi=0 



\j=0 

«2 + l 



1 «2+J / 

+ 1 Zj 'I 

7=0 V 



a 2 + 1 



- a 2 + zi + j - 1; v) - £(-a, + z x ; v)) 



(£(-ai - a 2 + j - 1; v) - f(-ai; v)) . 



fl2+lt^ V 7 

When v = this yields: 

ev ^=o ( ev a=o (^( -fl i» -«2)(Zi, z 2 ))) := ev^ (ev^ (f(-ai, -a 2 )(zi,z 2 ; 0))) 

(60) = _!_ + % J± . / 

a 2 + 1 ^) V ■/ /V «l+«2-J + 2 fli + l ( 

We next compute 



ev a=o ( ev zi=o -^2; v)(zi, Z 2 ))) 
5 



reg I reg 



+ev 



reg 
Z2=0 



a 2 - z 2 + 1 

[«2 - Z2]y-i 



(f(-ai -CL2+Z1+Z2-U v)-£(-ai +zi; v)) 



ev 



reg 
zi=0 



V 



2> 



Vi=l 



7 



- a 2 + zi + z 2 + 7 - 1; v) - £(-a, + z x ; v)) 



+ev 



reg 
Z2=0 



ev 



reg 



Z ^ 



[«2 - Z2];- 



(£(-ai - a 2 + zi + z 2 + 7 - 1; v) - £(-a, + zi; v)) 



ev reg . n | B ° , (£(-ai - a 2 + z 2 - 1; v) - £(-ai; v)) 



Z2=0 \a 2 + 1 

/ai+l 



+ev 



reg 



Z2=0 
«2+l 



Z B > 

Vi=i 



[a 2 - z 2 ]j- 



"2 + 1 / 

+ ^ B 7+a 1+ l 5 ;2 I 
7=1 V 



[a 2 - Z 2 ] y +ai 



(£(-ai - a 2 + z 2 + 7 - 1; v) - f (-ai ; v)) 
Res Z2=0 (£(-a 2 + z 2 + j; v)) 
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Q2 + 1 



— 7", t 



<a 2 + 1 



(£(-ai - a 2 + j - 1; v) - ^(-fli; v)) 



+(-l) ai+i a 1 !a 2 ! 



(ai + a 2 + 2)1 

since the only contribution to the residue comes from the term j = a\ + a 2 + 2. Since £(-a) 
£(-a; 0) = --f^f, this combined with (|5Q|) applied to = a 2 + 1 yields 



e\ 2 =o ( ev zi=o (^( -a i» -«2)(zi. z 2 ))) := ev^ (ev^ (f -«2! 0)(zi, z 2 ))) 

1 02+ 1 / 1 \ 



7=0 

+(-l) Ul+l ai!a 2 ! 



ai + a 2 — 7 + 2 ai + l a 2 + l 



(ai + a 2 + 2)! 
Combining ([6(1 and d6B yields 

^ ev (-ai,-a 2 ) := ev^"' 83 ™ (£(-ai, -a 2 )(zi,z 2 )) 

t „ /^2+ t\ c a 1 +a 2 -j+2 



(62) = -Vfi; 

«9 + 1 1 



«2 + 1 ^ \ J } a 1 +a 2 -j + 2 

+ I) Q,\ \a 2 ■ 



a.\ + 1 a 2 + 1 2 (ai + a? + 2) ! 

Renormalized multiple zeta values of depth 2 at non-positive arguments obtained this way, are 
rational linear combinations of Bernoulli numbers, and hence rational numbers. More generally, 
an inductive procedure on k carried out in the same spirit as the proof of Theorem 10 in 11551 shows 
that the renormalized multiple zeta values g E (sw' , Sk,v) are rational values at non-positive 
integer arguments s\, • ■ • , s k whenever v is rational. 

6.6.3. Multiple zeta values renormalized via Birkhoff decomposition. The renormalized multiple 
Hurwitz zeta values derived from a Birkhoff decomposition: 

<T Bir W..^; vi,...,v*) := y* Birk (o- Sun ®---®cr ShVk ) 

denoted by £ Birk (.si, . . . , s k , v) when vi = • • • = v k = v, satisfy stuffle relations 

£ Bkk (u*u';v) = £ B "\u;v){ B "\u';v). 

A striking holomorphy property arises at non-positive integer arguments [55] after implement- 
ing the diagonal map 5. 

Proposition 6.19. At non-positive integer arguments Si and for a rational parameter v, the map 

Z H> \l/*{0- SuV ®---®0- Sk , v )(z) 

defined in is holomorphic at zero. 
Consequently, 

^ Birk Oi, s k ; v) = lim ifr\<r suv ® • • • ® £r, foV ). 

z-*0 
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Let us compute double zeta values at non-positive integer arguments using Birkhoff decomposi- 
tion. Setting z\ = z 2 = z in (|59| ) leads to 



02 + 1 



Z[a 2 - z]j-i 
Bj (£(-ai - a 2 + 2z + j - 1; v) - £{-a y + z\ v)) 



J 



[«2 _ Z\a 2 +2 



+ - 



(«2 + 2) 



p £ |(/z + v) fll " z B^^Cy) (y + vT 2 dyj , 



Evaluating this expression at z = in a similar manner to the previous computation, yields: 



£ Birk (-a u -a 2 ;v) 



lim 

«2 + l 



'cf2 + l 
\j=0 



[02 ~ Z]j- 



(^(-ai - a 2 + 2z + j - 1; v) - £(-ai + z; v)) 



= £ 5, ^(" fl i - «2 + 7 - 1; v) - ^(-fli; v)) 



7=0 

+(-l) ai+1 ai!a 2 ! 



B 



ai+ci2+2 



2(a l +a 2 + 2)1 



When v = this yields [521: 



^k ( _ fli5 _ fl2) . = ^ ( _ au _ a2 . 0) 



1 



«2 + l 



+ 



"1 + - ;=0 



a 2 + 1\ 5 a 



H+ci2—j+2 



j j a l +a 2 - j + 2 



a\ + 1 a 2 + I 



+ (-l) Ql+1 ai!a 2 ! — — 



2(ai + a 2 + 2)! 



which coincides with (l62l) . 

Thus, renormalized double zeta values at non-positive integers obtained by two different meth- 
ods - using the symmetrized renormalized evaluator evQ en sym or a Birkhoff decomposition- coin- 
cide. 

Formula (l62l) yields the following table of values £(-a\, -a 2 ) for a\, a 2 e {0, . . . , 6} derived in 
H55J: 
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a-a, -b) 


a = 


a = 1 


a = 2 


a = 3 


a = 4 


a = 5 


a = 6 


b = 


3 
8 


i 

12 


i 

720 


i 

120 


n 

2 520 


i 

252 


i 

224 


b = 1 


1 

24 


1 

288 


1 

240 


19 


1 

504 


41 


1 

480 


10 080 


20160 


b = 2 


7 
720 


1 

240 


A 

u 


1 

504 


113 


1 

480 


307 


151200 


166 320 


b = 3 


1 

240 


1 

840 


1 

504 


1 


1 

480 


281 


1 

264 


28 800 


332 640 


b = 4 


11 

2 520 


1 

504 


113 


1 

480 


A 

u 


1 

264 


117 977 


151200 


75 675 600 


b = 5 


1 

504 


103 


1 

480 


1 

1232 


1 

264 


1 


691 


60480 


127 008 


65 520 


b = 6 


1 

224 


1 

480 


307 


1 

264 


117 977 


691 





166 320 


75 675 600 


65 520 



This table of values differs from the one derived in [|36l (see Table (|45l )) with which it however 
matches for arguments (a, b) with a + b odd and b ± and for diagonal arguments (-a, -a). 
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